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I.  Introduction 

Let  X be  a finite  set  of  real  points  and  let  f be  a function  defined  on  X 
(given  in  tabular  form)  which  we  desire  to  approximate.  Let  a = minfx:  x € X} 
and  b = max{x:  x € X}.and  for  any  function  g defined  on  X define  1 1 S 1 1-  x = max0s^x^l: 
x £ X}.  Finally,  let  TOL,  SMTH  and  N be  parameters  supplied  by  the  user  where 
TOL  is  a positive  number,  SMTH  is  a nonnegative  integer  and  N-l  is  an  integer  greater 
than  or  equal  to  SMTH.  In  this  setting,  our  algorithm  will  calculate  an  approximation, 
p,  to  f and  a set  of  points  (x.}._  C X with  a = x < x < ...  < x,  = b such  that  p 
restricted  to  [x.,  ] is  a polynomial  p^  G II  = {q:  q is  a real  algebraic 

polynomial  of  degree  <_  N-l},  p has  SMTH  continuous  derivatives  and  ||f-p||  ^ <_  TOL. 

In  what,  follows,  we  shall  use  the  notation  ||g||  g to  denote  max  {|g(t)|:  ,t  £ S} 
for  all  g defined  on  S and  S C X. 


II . The  Algorithm 

The  algorithm  begins  by  ctoosing  x^  to  be  the  largest  point  in  X such  that 

f\j 

(1)  [a,  x^]  H X contains  a\  least  max(2,  N + 1)  points,  and 

2 

(2)  If  p,  is  the  best  (unif*rm  or  L ) approximation  to  f from  IL.  , on 

1 N-l 

S^  = [a,  X ] 0 X then  ||i'  ” P q II  s — 

% ^ , . 

If  = b,  then  since  is  a oiecewise  polynomial  meeting  our  requirements, 

we  successfully  terminate  the  algorithm.  If  no  such  x^  exists,  the  algorithm 

fails  and  an  appropriate  error  message  is  generated.  Otherwise,  if  SMTH  = 0 

(i.e.,  we  only  require  the  approximation  to  be  continuous)  we  choose  the  right 

'V/ 

endpoint  of  the  first  subinterval,  x^,  to  be  x^.  If  SMTH  > 0,  in  order  to  add  to 
the  stability  of  the  algorithm  we  in  general  choose  x^  by  "backing  off"  from 

r\j 

x^  in  the  following  manner. 

We  first  examine  the  error  curve  f(x)  - p^(x)  and  find  those  points 

•••>  an  ^a»  X}-1  I"*  x at  which  relative  extrema  occur.  (Note  that  when 


interpolatory  constraints,  these  points  are  readily  available).  We  will  choose 

one  of  the  C^'s  to  be  x^.  The  motivation  for  choosing  x^  in  this  manner  is  that 

in  the  continuous  setting,  if  f is  differentiable  and  £ is  an  (interior)  relative 

extreme  point  of  f(x)  - p^(x)  then  f'(£)  - p^(£)  = 0 so  that  the  derivative  of 

p^  at  5 would  match  that  of  f at  £.  This  guarantees  that  when  we  smoothly  join 

the  next  piece  of  the  approximation  to  p^  at  £,  this  next  piece  will  closely  follow 

the  direction  of  f at  least  near  £.  If  we  merely  joined  the  second  piece  to  the 

first  at  x^  no  such  guarantee  can  be  made  and,  in  fact,  severe  oscillatory 

problems  tend  to  set  in.  Our  numerical  experience  indicates  that  the  procedure 

of  backing  off  from  x^  to  a smaller  x^  contributes  very  significantly  toward 

% 

the  stability  of  the  algorithm.  To  continue,  let  f’(£v)  be  the  derivative  of  the 

centered  quadratic  interpolation  of  f evaluated  at  £ . We  then  choose  x^  to  be 

% 

the  largest  £v  such  that  |f:’(£v)  - p^(Cv)|  < EPS,  where  EPS  is  a user-definable 

prescribed  tolerance,  or,  if  there  does  not  exist  such  a then  we  let  x^  be  the 

largest  Cv  at  which  | f ' ( ) - p^(£y)|  is  a minimum.  (In  our  implementation  of 

the  algorithm,  we  do  not  in  general  consider  all  of  the  relative  extreme  points  of 

f(x)  - p^(x)  in  (a,  x^]  0 X,  but  only  the  largest  N-SMTH-1  of  them.) 

We  continue  the  algorithm  by  finding  successive  intervals  [x, , x2],  [x2>  x3^> 

...»  [x  , b],  and  corresponding  polynomial  approximations  p , p , ...,  p g n 

m-i  z.  o m N-l 

to  f so  that  P^V!)=pij)(Vl)for  ^ = °»  •••»  SMTH  and  llf  - Pv  II  s £TOL, 

where  = [xy_^,  xy]  fj  X for  v = 2,  ...»  m,  xm  = b.  This  is  accomplished  as 

follows : 

Suppose  we  have  found  the  subintervals  [a,  x^,  [x1?  x2],  ...,  [x.  2 , x. 
and  the  corresponding  approximations  p^,  p2,  ...,  p^_  . Assume  further  that 
[Xj,  b]  C X contains  at  least  max(2,  N-SMTH)  points.  We  now  determine  an 
x^  and  a p^  meeting  the  above  requirements.  We  begin  by  choosing  x^  € X to  be 


the  largest  point  in  X which  satisfies 


4 


(1)  [x.  , x.]  H X contains  at  least  max(2,  N-SMTH)  points, 

(2)  If  p.  is  the  appropriate  best  approximation  to  f from  IT  , on 

i N-l 

S.  = [x.  , x.]  n X subject  to  the  constraint  that  p^(x.  ,)=p^(x.  , ), 

a l-l  i i l-l  i-I  l-l 

j = 0,1,...,  SMTH,  then  ||f  - p.  ||  < TOL. 

i bj 

<\j  % 

If  x.  = b,  we  set  x^  = x.  = b,  and  the  algorithm  is  successfully  terminated. 

If  no  such  X£  exists,  the  algorithm  fails  and  is  terminated.  Otherwise,  we 

choose  x.  completely  analogous  to  our  choice  of  x^  above. 

Finally,  we  consider  the  special  case  where  [x^_^,  b]  f)  X contains  fewer  than 

max(2,  N-SMTH)  points.  In  this  case  we  replace  x^_^  with  some  x^  ^ & X,  where 

x.  _ < x.  ,<x.  ,<b,  so  that  [x.  , , b]  D X contains  at  least  max(2,  N-SMTH) 
i-2  i-I  i-I  i-I  ’ 

points.  Specifically,  we  choose  to  be  a point  in  X closest  to  (b  - x^  2)/2 

which  satisfies 

(1)  [x.  ^ , b]  0 X contains  at  least  max(2,  N-SMTH)  points, 

(2)  If  p.  is  the  appropriate  best  approximation  to  f from  IT  , on 

1 N-l 

S = [x.  , b]  O X subject  to  the  constraint  that  pf^(x.  .)  = p.^  \x.  .), 

1 J.  1 1 — 1 1-1  1—1 


j = 0,  1,  . . . , SMTH,  then  ||f  ~ Pj.  II  s — T0L" 

Again,  if  we  can  find  such  an  x.  ^ then  the  algorithm  is  successfully 
terminated.  If  not,  the  algorithm  is  terminated  and  an  appropriate  error  message 
is  generated. 

N-l 

Remarks . Since  any  p.(x)  can  be  written  p^(x)  = £ a. (x  - x.  )^ , it  is  a simple 

j=0  1 1_1 

matter  to  meet  the  smoothness  constraints.  In  fact,  it  is  easy  to  generalize  this 
scheme  in  order  to  force  p^  to  satisfy  Hermite  interpolatory  constraints  at  several 
points.  Hence,  it  would  be  easy  to  modify  the  above  algorithm  so  that  we  could 


interpolate  the  function  being  approximated  and  its  derivatives  at  the  knots  if 
desired. 


5 


In  our  implementation  of  this  algorithm  we  have  modularized  the  approximation 

routine  so  that  we  may  use  a Remes-type  algorithm  for  computing  uniform  approximation 

subject  to  interpolating  constraints  [1],  and  we  use  Householder  transforms  to 

2 

compute  the  discrete  L approximations. 

Note  that  when  using  uniform  approximations  computed  by  the  Remes  algorithm 
it  is  in  general  not  necessary  to  compute  the  best  approximation  on  a particular 
subinterval  in  order  to  conclude  that  the  subinterval  is  too  long.  Indeed,  at 
each  iteration  of  the  Remes  algorithm  we  obtain  a lower  bound  for  the  error  of 
the  best  approximation  on  this  particular  subinterval.  Consequently,  if  during 
some  iteration  of  the  Remes  algorithm  this  lower  bound  is  greater  than  TOL,  we 
may  conclude  that  the  present  subinterval  is  too  large  and  terminate  the  Remes 
algorithm. 

. * % 

In  our  implementation  of  this  algorithm,  the  are  chosen  as  follows.  At 

r\j  m 

each  step  of  this  iterative  procedure  we  will  let  a be  the  current  largest  point 
in  X such  that  requirements  (1)  and  (2)  of  the  appropriate  algorithm  are  satisfied 
on  [x.^,  a]  n X,  and  we  will  let  % be  the  current  smallest  point  in  X such  that 
b > a and  requirement  (2)  of  the  appropriate  algorithm  fails  to  be  satisfied. 

We  initialize  this  process  by  computing  (or  attempting  to  compute)  the  appropriate 
best  approximation  on  [x^  b]  0 X.  If  this  approximation  satisfies  requirement 
(2),  then  we  set  = b and  we  are  done.  If  the  approximation  fails  to  satisfy 

*\j 

(2),  we  set  d = b.  Next,  let  t = min{x  £ X:  [x.  x]  f)  X contains  at  least 
max(2,  N-SMTH)  points}.  If  the  approximation  on  [x.._^,  t]  P)  X fails  to  satisfy 
(2)  then  the  algorithm  cannot  meet  the  desired  accuracy  and  fails.  Otherwise, 
we  set  a = t. 

In  general,  we  proceed  as  follows.  We  let  t = inf{x  £ X:  (i>  - a)/2  <_  x < B}. 

If  this  set  is  empty,  we  set  t = sup{x  6 X:  a £ x <_  (fe  - a)/2}.  If  t = a then 

*\j  % , 

this  procedure  has  converged  and  we  set  x^  = t = a.  Otherwise,  we  compute  (or 
attempt  to  compute)  the  appropriate  approximation  on  [x._^,  t]  D X.  If  this 


approximation  satisfies  (2)  then  we  set  a = t.  If  this  approximation  fails  to 


satisfy  (2)  then  we  set  b ~ t.  We  continue  this  process  until  b - a is  less 


than  some  user  definable  prescribed  tolerance,  at  which  point  we  accept  a as  a 


and  terminate  this  procedure.  We  compute  the  x.  in  a 


good  approximation  to  x 


manner  analogous  to  the  above 


In  an  attempt  to  accelerate  the  convergence  of  the  above  scheme  when  using 


the  best  uniform  approximation  operator,  we  have  tried  to  take  advantage  of  the 


call  it  SMLERR,  and  corresponding  to  each  b we  know  a lower  bound 


call  it  BIGERR 


above  scheme  by  first  setting  a = (BIGERR  - TOL)/(BIGERR  - SMLERR)  and  then 


setting  t = infix  £ X:  aa  + (1  - a)b  < x < b}  or,  if  this  set  is  empty,  we  set 


Hence,  if  SMLERR  is  very  close  to  TOL,  t will  be  chosen  close  to  a.  Our  numerical 


experience  has  shown  that  this  procedure  only  works  well  when  approximating 


uniformly  smooth  functions  and  that  this  algorithm  cannot  be  significantly 


improved  by  allowing  the  Remes  algorithm  to  run  to  completion  in  order  to  obtain 


the  true  error  of  approximation  for  BIGERR 


III.  Numerical  Results 


This  algorithm  has  been  implemented  as  FORTRAN  programs  and  has  been  tested 


on  Colorado  State  University's  CDC  CYBER  1?2.  In  the  appendix  we  give  a listing 


of  the  algorithm  using  the  best  uniform  approximation  operator.  By  using  a least 


squares  routine  in  place  of  the  REMES  subroutine  (and  those  subroutines  called 


only  by  REMES:  SOLVE,  DIVIDF , ZEROFD,  BPOLY,  TRANS)  the  least  squares 


version  of  this  algorithm  is  readily  obtained.  Indeed,  we  essentially  replaced 


these  six  subroutines  by  four  subroutines:  ASET,  HOUSE,  TOLCHK  and  FIX.  ASET 


sets  up  the  overdetermined  system  to  be  solved  in  the  L 


4 


7 

amount  of  care  needs  to  be  taken  to  insure  that  the  approximation  on  each  succeeding 
interval  has  the  desired  smoothness  at  the  common  endpoints.  Thus,  if  we  are 
currently  considering  the  interval  [x^_.p  x,.]  and  Pi-i  aPProx:'-ma't i°n 

that  was  accepted  for  [x. x.  j]  then  we  form  the  following  system: 


N-l 

I c 


(x-x.  )V=f (x)-  l 
L 1-1  u „ 


k P^x.^)  fx£  [x^.x.]  n X if  k = -1 

- — l.  ; (x-x.  ) for  < 

v=k+l  V 1 1 v=0  V-  1 [x£  (x^i.x.DD  X if  k > -1 

-1 

where  k = SMTH  = smoothness  desired  at  x.  . and  £ is  to  be  replaced  by  0 

1-1  v=0 

whenever  it  occurs.  HOUSE  is  any  good  over determined  least  squares  solver. 

We  use  one  based  on  Householder  transforms  (a  simple  version  since  our  system 
will  always  be  nondegenerate).  TOLCHK  compares  the  maximum  absolute  value  of 
f minus  the  approximation  returned  by  HOUSE  with  TOL  and  also  finds  the  set  of 
N - k - 1 last  "extreme"  values  in  [x^_^,  x^] . That  is,  points  at  which 
| f ( x ) - p ( x ) | is  a local  maximum  relative  to  points  on  each  side  of  x.  This 

is  for  the  backing  off  from  x.  procedure  to  increase  stability.  FIX  simply 

1 N v-1 

converts  the  coefficients  {c  } of  a polynomial  of  the  form  / c (x-x.  , ) 

v r , v l-l 

v=l 


N 


into  coefficients  {b  } for  the  same  polynomial  written  in  the  form  £ b xV  1. 

v=l  v 

Ixl 

As  examples,  using  the  best  uniform  operator  the  functions  e 1 on  [-1,1], 
|sin(x)|  on  [— u,ir]  qnd  /x  on  [0,2]  were  approximated  on  200  equally  spaced 
points  with  TOL  = .0,1,  SMTH  = 2,  and  N = 6.  Each  of  these  examples  is  difficult 
to  approximate  by  polynomials  near  x = 0.  Consequently,  the  algorithm's  ability 
to  automatically  decrease  the  length  of  the  subintervals  near  x = 0 and  then 
recover  by  lengthening  them  for  x > 0 is  tested.  Below  is  a table  listing  knot 
locations  (subinterval  endpoints)  and  the  CPU  time  in  seconds  used  to  compute 
the  piecewise  polynomial  approximation. 





Because  uniform  approximations  weight  each  data  point  equally,  they  are 
particularly  suited  for  approximating  precise  mathematical  functions  or  for 
approximating  data  with  noire  levels  that  are  very  small  relative  to  the 
desired  accuracy.  However,  when  these  algorithms  are  used  to  approximate  functions 
containing  considerable  levels  of  noise,  in  general  the  approximations  tend  to 
follow  the  noise  patterns  more  than  is  desirable  so  that  near  abrupt  changes  in 
the  data,  the  approximations  tend  to  begin  to  oscillate  and  generally  it  requires 
several  subintervals  to  dampen  these  oscillations.  (Some  counterexamples  from 
engineering  applications  have  been  found  which  show  that  this  is  not  always  the 
case . ) 

For  example,  we  were  given  some  experimental  data  involving  the  release  of 
bitumen  from  oil  shale  heated  to  a constant  temperature  as  a function  of  time. 
Because  relatively  few  data  points  were  available,  we  filled  in  the  gaps  between 
the  data  points  by  linear  interpolation  so  that  we  approximated  on  200  equally 
spaced  points.  Figure  1 is  a plot  of  the  data  being  approximated  and  the  approxi- 
mation generated  using  uniform  approximation  with  N = 6,  SMTH  = 2,  and  TOL  = 2.5. 


_ 


- 


For  some  applications  it  might  be  preferable  to  be  able  to  choose  ahead 
of  time  a particular  knot  location  and  to  specify  interpolatory  constraints  at 
these  knots.  For  example,  when  approximating  [sin(x)j  it  might  be  advantageous 
to  force  a knot  to  be  located  at  x = 0 and  to  require  p(0)  = 0,  p'(0)  = -1  as 
x approaches  0 from  the  left,  and  p'(0)  = tl  as  x approaches  0 from  the  right. 

The  modifications  to  the  algorithm  as  given  necessary  to  accomplish  this  would 
not  be  extensive  or  difficult. 

It  should  also  be  remarked  that  this  algorithm  can  be  used  with  linear 
approximating  families  other  than  polynomials.  For  example,  any  generalized 
Haar  subspace  of  high  enough  order  (to  allow  Hermite  interpolating  constraints) 
can  be  used  with  the  best  uniform  approximation  operator.  In  addition,  one 
could  force  not  only  knots  to  occur  at  various  prescribed  points  but  could  also 
vary  the  smoothness  at  these  knots  as  well  as  change  the  form  of  the  approximating 
families  at  these  -knots. 

Recently,  John  Rice  has  described  an  adaptive  piecewise  polynomial  algorithm 
which  uses  local  Hermite  interpolation  instead  of  uniform  approximation  on  each 
subinterval  to  obtain  each  polynomial  piece  [2].  His  algorithm  requires  that  the 
first  SMTH  derivatives  or  approximate  derivatives  of  f be  available  in  order  to 
compute  approximations  which  have  SMTH  continuous  derivatives.  Rice's  adaptive 
strategy  also  differs  from  ours;  he  uses  a bisection  strategy  which  can  be 
described  as  follows.  First,  compute  the  Hermite  interpolating  polynomial  to  f 
on  [a,  b].  That  is,  compute  the  polynomial  which  interpolates  f and  its  first 
SMTH  derivatives  at  both  endpoints,  as  well  as  interpolating  f at  k evenly 
distributed  points  in  (a,  b),  where  k = DEGREE  - 2* SMTH  - 1,  and  DEGREE  is  the  degree 
of  the  approximating  polynomial.  Next,  measure  the  error  of  approximation  using 
any  preselected  1/  norm  (p  >_  1).  If  the  error  is  less  that  TOL,  the  algorithm 
terminates,  otherwise  bisect  [a,  b]  into  two  subintervals  and  check  to  6ee  if 
the  Hermite  interpolating  polynomial  on  [a,  (a  + b)/2]  differs  from  f in  norm  by 


! 

I 


t 


•I  • 
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no  more  than  TOL.  Meanwhile,  place  the  subinterval  [(a  + b)/2,  b]  in  a "stack" 
to  be  processed  later.  Continue  bisecting  subintervals,  placing  the  right  half 
of  the  interval  on  the  top  of  the  stack  to  be  processed  later  and  check  the  left 
half  to  see  if  the  error  of  approximation  by  the  Hermite  interpolating  polynomial 
is  less  than  TOL.  When  a subinterval  and  its  associated  approximations  are  found 

. 

which  meet  the  desired  tolerance,  we  accept  this  approximation  as  a "piece"  of  the 
piecewise  polynomial  approx imation  and  continue  the  algorithm  by  removing  the 
subinterval  which  is  currently  at  the  top  of  the  stack  and  repeating  the  above 
on  it  unless  the  stack  is  empty,  at  which  point  terminate  the  algorithm. 

Rice's  routine  requires  the  value  of  the  function  and  its  derivatives  at 
very  many  points  in  [a,  b],  even  though  it  may  not  actually  access  these  values. 

His  algorithm  is  particularly  suited,  then,  for  approximating  precise  mathematical 
functions  for  which  this  data  is  readily  available.  Our  algorithm  seems  more 
suitable  for  obtaining  approximations  of  functions  given  in  tabular  form  (as  in 
many  engineering  applications)  since  our  algorithm  does  not  require  any  information 
about  derivatives,  and  no  special  configuation  of  the  data  is  needed  for  our 
adaptive  strategy  to  work.  See  the  forthcoming  paper  by  Andrews,  Hull  and  Taylor 
for  further  comments  on  suitability  of  this  algorithm  for  this  type  of  application. 

By  using  interpolatory  constraints  on  both  endpoints  of  each  subinterval, 
and  interpolating  the  values  of  the  derivative  of  f at  these  points,  we  could 
modify  our  algorithm  to  use  a bisection  adaptive  strategy.  However,  in  order  to 
take  greatest  advantage  of  our  use  of  the  best  approximation  operators  as  opposed 
to  the  Hermite  interpolation  operator,  we  felt  that  it  was  preferable  to  use  the 
majority  of  our  degrees  of  freedom  inside  subintervals  instead  of  using  them  up 
on  interpolation  requirements  at  the  endpoints.  Essentially,  in  our  algorithm  we 
traded  a somewhat  faster  run  time  (since  computing  best  approximations  is  very  slow 
relative  to  computing  interpolatory  polynomials)  for  generally  fewer  knots  and  the 
freedom  from  needing  to  know  derivative  information. 
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APPENDIX 

Here  we  give  listings  for  the  two  algorithms  described  in  this  paper  with  the 
uniform  approximation  listing  first  and  the  least  squares  approximation  listing 
second.  In  both  cases  we  include  a driver  and  sample  output  where  the  output 
corresponds  to  figure  1 and  figure  2,  respectively.  Also,  we  list  the  input 
of  these  two  runs  between  the  listing  of  the  code  and  the  output. 
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UNIFORM  ADAPTIVE  CUhVt  FITTING  PROGRAM  t SAMPLE  RUN 


INPUT  1 

6 2 

2.50 

( N,  SMTH,  TOL  ) 

3.0 

0.0 

( XTABLE , FTABLE) 

5.0 

1.3 

7.0 

3.4 

11.0 

5.2 

13.0 

6.0 

15.0 

14.4 

17.5 

21 .4 

20.0 

27.4 

22.5 

50.9 

25.0 

49.3 

27.5 

47.5 

30. 0 

51.5 

35.0 

36.5 

40.0 

27.9 

50.0 

9.4 

60.0 

4.2 

OUTPUT  i 

INTERVAL  NUMBER  1 WHICH  BEGINS  AT  .3000000000000000E+01 
AND  ENDS  AT  .219045226l3065l8E+02  CONTAINS  67  POINTS. 

THE  COEFFICIENTS  OF  LEST  APPROXIMATION  IN  THIS  INTERVAL  ARE 

c(  1)  = -.1634025000676115E+02 

C(  2)  = .1181130531149097E+02 

C(  3)  = - . 263714627291 8883E+ 01 
C(  4)  = .2593743838385461E+00 

C(  5)  = -.1129567576998175E-01 
C(  6)  = .1865450590541398E-O3 

THE  ERROR  OF  APPROXIMATION  IN  THIS  INTERVAL  IS  . 1969452589748130E+01 . 


INTERVAL  NUMBER  2 WHICH  BEGINS  AT  .2190452261306518E+02 
AND  ENDS  AT  . 270603015075374 7E+02  CONTAINS  19  POINTS. 

THE  COEFFICIENTS  OF  BEST  APPROXIMATION  IN  THIS  INTERVAL  ARE 

C(  1)  - .7960326211065352E+06 

C(  2)  *=  -.1656777306495747E+06 
C(  3)  » .1375385536550829E+05 

C(  4)  * 56927784 74236584E+03 

C(  5)  = . 11748807 89994893E+ 02 

C ( 6)  = - .96 7283885240641 8E -01 

THE  ERROR  OF  APPROXIMATION  IN  THIS  INTERVAL  IS  . 22948744593741 14E+01 . 


INTERVAL  NUMBER  3 WHICH  BEGINS  AT  .27060301 5075374 7E+ 02 
AND  ENDS  AT  .2849246231155757E+02  CONTAINS  6 POINTS. 

THE  COEFFICIENTS  OF  BEST  APPROXIMATION  IN  THIS  INTERVAL  ARE 


l 


29 


C(  1)  = - .^5555223084l8542E+07 
C(  2)  = . 781533098^931 737E+ 06 

C(  3)  = -.5356190391737269E+05 
C(  4)  = . 1 83305191 31 2759^E+C4 

c(  5)  = -.3132567622722183E+02 

C(  6)  = .2138556370651710E+00 

THE  ERROR  OF  APPROXIMATION  IN  THIS  INTERVAL  IS  . 2350077767167932E+01 . 


INTERVAL  NUMBER  4 WHICH  BEGINS  AT  . 2849246231155757E+02 
AND  ENDS  AT  . 3078391959798978E+02  CONTAINS  9 POINTS. 

THE  COEFFICIENTS  OF  BEST  APPROXIMATION  IN  THIS  INTERVAL  ARE 

c(  1)  = .1040230085838158E+07 

C(  2)  = -.1649042600538041E+06 
C(  3)  = . 1043144 021 371 787F+ 05 

C ( 4)  = -.3291270419910106E+03 
C(  5)  = .5179689635444362E+01 

C(  6)  = 325294 2438648798E-01 

THE  ERROR  OF  APPROXIMATION  IN  THIS  INTERVAL  IS  .2350077?448l6l90E+01 . 


INTERVAL  NUMBER  5 WHICH  BEGINS  AT  . 3078391959798978E+02 
AND  ENDS  AT  .3794472361809017E+02  CONTAINS  26  POINTS. 

THE  COEFFICIENTS  OF  BEST  APPROXIMATION  IN  THIS  INTERVAL  ARE 

c(  1)  = -.232982069020337^+06 

C(  2)  = .3289571^88642064E+05 

C(  3)  = -.1849503807458532E+04 
C(  4)  = .5l77682249693^07E+02 

C(  5)  = -.7218440333123937E+00 
C(  6)  = .4009674821594894E-02 

THE  ERROR  OF  APPROXIMATION  IN  THIS  INTERVAL  IS  . 2347864915846003E+01 . 


INTERVAL  NUMBER  6 WHICH  BEGINS  AT  .379447236l8090l?E+02 
AND  ENDS  AT  .4310050251256257E+02  CONTAINS  19  POINTS. 

THE  COEFFICIENTS  OF  BEST  APPROXIMATION  IN  THIS  INTERVAL  ARE 

C(  1)  = -,6827^3^677450825E+05 
C(  2)  = .7613513^11 31 0967E+ 04 

c(  3)  = -.3373191291063522E+03 
C(  4)  = .7428809770578283E+01 

c(  5)  ■ -.8135813032576H1E-01 
C(  6)  * .35l*5^53066058774E-03 

THE  ERhOR  OF  APFhOXI NATION  IN  THIS  INTERVAL  IS  . 24l602062504944lE+0l . 


INTERVAL  NUMBER  7 WHICH  BEGINS  AT  .431 0050251256257E+02 
AND  ENDS  AT  . 5999999999999955E+02  CONTAINS  60  POINTS. 

THE  COEFFICIENTS  OF  BEST  APPROXIMATION  IN  THIS  INTERVAL  ARE 

L 


.1 33361 6589900339E+ 05 
-.1192965606757301E+04 
. 419387066 7632291E+02 
-.7183879075542343E+00 
.593^7961 53467731E-02 
-.1861388516033458E-04 


THE  ERROR  OF  APPROXIMATION 


THIS  INTERVAL  IS  . 24l60206248748l8E+01 . 
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LEAST  SQUARES  ADAPTIVE  CURVE  FITTING  PROGRAM  i SAMPLE  RUN 


INPUT  i 


6 2 

3.00 

( N,  SMTH,  TOL  ) 

3.0 

0.0 

(XTABLE,  FTABLE) 

5.0 

1.3 

7.0 

3.4 

11.0 

5.2 

13.0 

6.0 

15.0 

14.4 

17.5 

21.4 

20.0 

27.4 

22.5 

50.9 

25.0 

49.3 

27.5 

47.5 

30.0 

51.5 

35.0 

36.5 

40.0 

27.9 

50.0 

9.4 

60.0 

4.2 

OUTPUT  1 

INTERVAL  NUMBER 

1 WHICH 

BEGINS  AT  .3000000000000000E+01 

AND  ENDS  AT  . 2l857l42857l4278E+02  CONTAINS  133  POINTS. 

THE  COEFFICIENTS  OF  BEST  APPROXIMATION  IN  THIS  INTERVAL  ARE 

C(  1)  = -.3189942743572442E+02 
C(  2)  = .1812103641220483E+02 

c(  3)  = -.3560852823483785E+01 
C(  4)  = .3242284790553871E+00 

c(  5)  = 135023490206 85 14E-01 

C(  6)  = .2158257913951744E-03 

THE  ERROR  OF  APPROXIMATION  IN  THIS  INTERVAL  IS  . 287469301288797E+01 . 


INTERVAL  NUMBER  2 WHICH  BEGINS  AT  . 2185714285714278E+01 
AND  ENDS  AT  . 252857l42857l4llE+02  CONTAINS  25  POINTS. 

THE  COEFFICIENTS  OF  BEST  APPROXIMATION  IN  THIS  INTERVAL  ARE 

C(  1)  = .7525331567081399E+06 

C(  2)  = -.1568254831136037E+06 
C(  3)  - .1303544628337352E+O5 

C(  4)  * -.54021 744 34562294E+03 
C(  5)  = .H16296306156516E+02 

C(  6)  = -.9202040460845584E-01 

THE  ERHOR  OF  APPROXIMATION  IN  THIS  INTERVAL  IS  . 2824472463024449E+01 . 


INTERVAL  NUMBER  3 WHICH  BEGINS  AT  . 2528571428571411E+02 
AND  ENDS  AT  .3557142857142844E+0 2 CONTAINS  73  POINTS. 

THE  COEFFICIENTS  OF  BEST  APPROXIMATION  IN  THIS  INTERVAL  ARE 


LI 


45 


C(  1) 
C(  2) 


c(  5) 

C(  6) 


.6867241093987040E+05 
1 1211 369880750 14E+ 05 
. 7277741 519933406E+03 
-.2346667791438097E+02 
• 3759555635732301E+00 
-.2395099334129056E-02 


THE  ERROR  OF  APPROXIMATION  IN  THIS  INTERVAL  IS  . 2985086?50638296E+01 . 

INTERVAL  NUMBER  4 WHICH  BEGINS  AT  . 3557142857142844E+02 
AND  ENDS  AT  . 5399999999999977E+02  CONTAINS  130  POINTS. 

THE  COEFFICIENTS  OF  BEST  APPROXIMATION  IN  THIS  INTERVAL  ARE 

C(  1)  = 844662809591 8451E+04 

C(  2)  = . 9702467401 092472E+ 03 

C(  3)  = -.4389461668285185E+02 
C(  4)  = ,9842698l9824i457E+00 

C(  5)  = -.1096799498228268E-01 
C(  6)  = . 48633^040126 7810E-04 

THE  ERROR  OF  APPROXIMATION  IN  THIS  INTERVAL  IS  . 14647851 27359050E+ 01 . 

INTERVAL  NUMBER  5 WHICH  BEGINS  AT  . 5399999999999977E+02 
AND  ENDS  AT  . 5999999999999977E+02  CONTAINS  43  POINTS. 

THE  COEFFICIENTS  OF  BEST  APPROXIMATION  IN  THIS  INTERVAL  ARE 

C(  1)  o -.1771038129739963E+07 
C(  2)  = .1591827758229449E+06 

C(  3)  = -.5718444554639893E+04 
C(  4)  = .1026330177930108E+03 

C(  5)  = -.9202908099941851E+00 
C(  6)  = .3298242672003487E-02 

THE  ERROR  OF  APPROXIMATION  IN  THIS  INTERVAL  IS  .9836415372650436E+00 . 


